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We discuss the non-Gaussian contribution to the power spectrum covariance of cosmic microwave 
background (CMB) anisotropics resulting through weak gravitational lensing angular deflections and 
the correlation of deflections with secondary sources of temperature fluctuations generated by the 
large scale structure, such as the integrated Sachs- Wolfe effect and the Sunyave-Zel'dovich effect. 
This additional contribution to the covariance of binned angular power spectrum, beyond the well 
known cosmic variance and any associated instrumental noise, results from a trispectrum, or a four 
point correlation function, in temperature anisotropy data. With substantially wide bins in multipole 
space, the resulting non-Gaussian contribution from lensing to the binned power spectrum variance is 
insignificant out to multipoles of a few thousand and is not likely to affect the cosmological parameter 
estimation with acoustic peaks and the damping tail. The non-Gaussian contribution to covariance, 
however, should be considered when interpreting binned CMB power spectrum measurements at 
multipoles of a few thousand corresponding to angular scales of few arcminutes and less. 



I. INTRODUCTION 

The applications of cosmic microwave background 
(CMB) anisotropy measurements are well known Qj; its 
ability to constrain most, or certain combinations of, pa- 
rameters that define the currently favorable cold dark 
matter cosmologies with a cosmological constant has 
driven a wide number of experiments from ground and 
space. The advent of high sensitivity CMB anisotropy 
experiments with increasing capabilities to detect fluc- 
tuations over a wide range of scales now suggest the 
possibility that anisotropy power spectrum at small an- 
gular scales will soon be measured. At angular scales 
corresponding to few arcminutes and below, fluctuations 
are mostly dominated by secondary effects due to local 
large scale structure (LSS) between us and the recombi- 
nation. Additionally, important non-linear second order 
effects such as the weak gravitational lensing of CMB 
leaves important imprints that can in return be used as 
a probe of cosmology or astrophysics related to evolution 
and growth of structures. 

The increase in sensitivity of current and upcoming 
CMB experiments also suggest the possibility that non- 
Gaussian signals in the CMB temperature fluctuations 
may be detected and studied in detail. The devia- 
tions from Gaussianity in CMB temperature fluctuations 
arise through two scenarios: the existence of a primor- 
dial non-Gaussianity associated with initial fluctuations 
and the creation of non-Gaussian signals through non- 
linear mode-coupling effects related to secondary contri- 
butions. In currently favored cosmologies with adiabatic 
initial conditions, the primordial non-Gaussianity is non- 
existent or insignificant 0. This leaves the secondary 
contributions, such as the Sunyave-Zel'dovich effect |3|] 
due to inverse- Compton scattering of CMB photons via 
hot electrons, as the main source of non-Gaussianity. The 
existence of non-Gaussian fluctuations in temperature 
can be directly measured through higher order correla- 



tions, such as a three-point function or a bispectrum in 
Fourier space The detection of non-Gaussianities at 
the three-point level can be optimized through the use of 
special statistics and matched filters |(| and through cer- 
tain physical aspects associated with secondary effects, 
such as frequency dependence J7|. 

An additional effect due to non-Gaussianity include a 
contribution to the four-point correlation function, or a 
trispectrum in Fourier space, of CMB temperature fluc- 
tuations The four point correlations are of spe- 
cial interest since they quantify the sample variance and 
covariance of two point correlation or power spectrum 
measurements |1G,0]. Thus, to properly understand the 
statistical measurements of CMB anisotropy fluctuations 
at the two point level, a proper understanding of the 
four point contributions is needed. Similarly, studies of 
the ability of CMB power spectrum measurements to 
constrain cosmology have been based on a Gaussian ap- 
proximation to the sample variance and the assumption 
that covariance is negligible. If there are significant non- 
Gaussian contributions from the four point level that con- 
tribute to the power spectrum covariance, then it could 
affect the conversion of power spectrum measurements 
to estimates on cosmological parameters. Given the high 
precision level of cosmological parameter measurements 
expected from CMB, a careful consideration must be at- 
tached to understanding the presence of non-Gaussian 
signals at the four point level. Thus, the basic goal of 
this paper is to understand to what extent Gaussian as- 
sumption on CMB power spectrum covariance remains 
to hold when non-Gaussian contributions are included. 

As discussed in previous studies (e.g., jl), one of the 
most important non-linear contribution to CMB temper- 
ature fluctuations is weak gravitational lensing. Simi- 
lar to the weak lensing contribution to CMB anisotropy 
at the three-point level, the contributions to the four 
point level results from the non-linear mode-coupling na- 
ture of weak lensing effect and the correlation between 
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weak lensing angular deflections and secondary effects 
that trace the same large scale structure. The trispec- 
trum due to lensing alone is studied in ||] and the same 
trispectrum, under an all-sky formulation, is also consid- 
ered in ||. Here, we focus on the contribution of the 
trispectrum to the power spectrum covariance which was 
not considered in previous works. We also include the 
trispectrum resulting from the correlation between lens- 
ing and secondary effects such as the interagted Sachs- 
Wolfe (ISW; ||) effect and the thermal SZ § effect. 
The latter effect has now been imaged towards massive 
galaxy clusters where the temperature of the scattering 
medium can be as high as 10 keV producing temperature 
fluctuations of order 1 mK |l3|] . 

In general, we do not expect secondary contributions 
such as the thermal SZ effect to be an important non- 
Gaussian contribution to temperature fluctuations, since 
its signal can be easily separated from the thermal CMB 
spectrum based on multifrequency information . Still, 
we consider its correlation with lensing as a source of co- 
variance as certain experiments, especially at small angu- 
lar scales, may not have the adequate frequency coverage 
for a proper separation. In a previous paper, we discussed 
the non-Gaussian covariance resulting from SZ alone is 
discussed in [Q, where we also considered the effect of 
full-covariance, compared to Gaussian variance assump- 
tion, on the estimation of parameters related to the SZ 
effect. As discussed there, due to highly non-Gaussian 
behavior of the SZ signal resulting from its dependence 
on massive halos such as galaxy clusters, the determi- 
nation of parameters that define the SZ contribution is 
affected by the presence of non-Gaussian contribution to 
the covariance. 

In §||, we introduce the basic ingredients for the 
present calculation. The CMB anisotropy trispectra due 
to weak lensing and correlations between weak lensing 
and secondary effects are derived in §|lj]. In §|y|, we cal- 
culate the CMB power spectrum covariance due to weak 
lensing and discuss our results in §|v|. In § fv| we conclude 
with a summary. 



II. GENERAL DERIVATION 

Large-scale structure between us and the last scatter- 
ing surface deflects CMB photons propagating towards 
us. Since lensing effect on CMB is essentially a distri- 
bution of photons, from large scales to small scales, the 
resulting effect appears only in the second order [g. In 
weak gravitational lensing, the deflection angle on the 
sky given by the angular gradient of the lensing poten- 
tial, 5(h) = V0(n), which is itself a projection of the 
gravitational potential, $ (see e.g. ||i~5|| ), 

<Km) = -2 f° dr dAi \°~ ( r) M r, rhr) . (1) 
J d A (r)d A {r ) 

The quantities here are the conformal distance or look- 
back time, from the observer, given by 
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FIG. 1. Power spectrum of the CMB temperature 
anisotropies in the fiducial ACDM model. We have also 
shown the weak lensing, integrated Sachs- Wolfe and the SZ 
secondary contributions to angular power spectrum. 
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and the analogous angular diameter distance 



d A = H^fl^ 2 smh(H Q^r) 



V2„ 



(2) 



(3) 



with the expansion rate for adiabatic CDM cosmological 
models with a cosmological constant given by 

h 2 = Hi [n m (i + zf + n x (i + z) 2 + n A ] • (4) 

Here, Hq can be written as the inverse Hubble distance 
today Hq 1 = 2997.9ft~ 1 Mpc. We follow the conventions 
that in units of the critical density 3Hq/8itG, the con- 
tribution of each component is denoted i — c for the 
CDM, b for the baryons, A for the cosmological constant. 
We also define the auxiliary quantities fi TO = Q c + Qb 
and fix = 1 — J2i which represent the matter density 
and the contribution of spatial curvature to the expan- 
sion rate respectively. Note that as VLk — * 0, d A — > r and 
we define r(z — oo) = tq. Though we present a general 
derivation of the trispectrum contribution to the covari- 
ance, we show results for the currently favorable ACDM 
cosmology with Of, = 0.05, O m = 0.35, ft a = 0.65 and 
h = 0.65. 

The lensing potential in equation [l] can be related to 
the well known convergence generally encountered in con- 
ventional lensing studies involving galaxy shear p5[ 
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1 , 
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(rh) 



(5) 



r ° dr d A (r)d A (r 
d A {r Q ) 



V^$(r, rhr) 



where note that the 2D Laplacian operating on $ is a 
spatial and not an angular Laplacian. Though the two 
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terms k and (j> contain differences with respect to radial 
and wavenumber weights, these differences cancel with 
the Limber approximation [^6|. The spherical harmonic 
moments of these two quantities are simply proportional 
with multiplicative factors in I 



1(1 + 1) 



J dny," l *(n)0(n) 



where 



W len (k,r) 



d 3 U 

^s(k)Yr(k)i i r(k), 



drW Xcn {k ir )ji{kr) , 



(6) 



p( s dA{r Q - r) 
k J {T> d A (r)d A (r ) 



(7) 



Here, we have used the Rayleigh expansion of a plane 
wave 



(8) 
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and the fact that V 2 V; m 



-1(1 + l)Y, m . 

Note that the cosmological Poisson equation relates 
fluctuations in the density field, given by 8 in equation 
to the fluctuations in the gravitational potential, <i>: 



i + 3§n K 



G(r) 



8(k,0), (9) 



with the growth function, S(k, r) — G(r)8(k,0) |17|, given 
by 



G(r) cx 



H(r) 



Jz(r) 



dz'{l + z') 



Hn 



H(z') 



(10) 



Note that in the matter dominated epoch G cx a = (1 + 
Expanding the lensing potential to multiple moments, 

= £> m iT(A). (11) 

we can write its correlation function 

(<Hn)0(m)) = E Wimh'm') Y r • (12) 

The power spectrum of lensing potentials is now given 
through 



Cf 



k 2 dkP{k)l\ CU (k)l\ en (k) 



(13) 



(14) 



Here, we have introduced the power spectrum of den- 
sity fluctuations 



(<S(k)<J(k')) = (27r)-M u (k + k')P(fc) 



where 8 is the Dirac delta function and 



k 3 P(k) 
2tt 2 



n+3 



(15) 



(16) 



in linear perturbation theory. We use the fitting formulae 
of in evaluating the transfer function T(k) for CDM 
models. Here, Sh is the amplitude of present-day density 
fluctuations at the Hubble scale; with n — 1, we adopt 
the COBE normalization for 8 H M of 4.2 x 10" 5 , consis- 
tent with galaxy cluster abundance [BOf, with erg = 0.86. 
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FIG. 2. Power spectra (a) of lensing angular deflections 
and (b) lensing-isw and lensing-sz cross-correlation. In (a), 
we show the lensing deflection power spectrum under linear 
perturbation theory description of the matter fluctuations and 
using the halo model. As shown, the lensing deflection power 
peaks at multipoles of ~ 40 and act as an effective window 
function which smooths the CMB power spectrum. In (b) , the 
lensing-SZ correlation is calculated with the halo approach 
while lensing-ISW correlation follows the use of linear theory 
dark matter power specttrum. 
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Note that an expression of the type in equation fll4| ) 
can be evaluated efficiently with a version of the Limber 
approximation fl6| , called the weak coupling approxima- 
tion Ml, such that 



drW x {k,r)ji{kr) 



W x (kJ/k) f drji(kr) (k^- > 1) 
Jo 1 W 

RY[{l + l)/2] 



effect, due to the line of sight motion of electrons, is sec- 
ond order in density fluctuations, there is no first order 
cross-correlation with the lensing potentials. The ISW 
|l2| effect results from the late time decay of gravita- 
tional potential fluctuations. The resulting temperature 
fluctuations in the CMB can be written as 



W x (k,l/k) 



2k T[(7 + 2)/2] 



(17) 



For I ^> 1, the ratio of gamma functions goes to 
simplifying further. Using a change of variables such that 
cLa = l/k, we obtain an approximation for the power 
spectrum of lensing potentials as 



Cf 



k 2 dkp{k)i\ cn {k)ir{k) 



(18) 



dr 



Since the same large scale structure responsible for de- 
flections in CMB photons produce contributions to the 
anisotropies through other effects, there is a correlation 
between the deflection potential and secondary sources 
of temperature fluctuations. As in the power spectrum 
of delection potential, using statistical isotropy, we write 
the correlation as 



(<^L«L> = bf - 



-cT K 



- / k 2 dkp{k)if{k)ir{k), 



(19) 



Here. 



dr 
"a 



T s (n) 



dA ' 



W 



lcn 



dA / \dA 



(20) 



and we have used equation (^|) to relate the power spee- 
dy and the K-secondary cross power 



trum bf of Refs. 

spectrum of Ref . p2| . The last line represents the Limber 
approximation and we have assumed that the secondary 
anisotropies are related to the first order fluctuations in 
the density field projected along the line of sight, 



Hm 



d 3 k 



^(k)7f(fc)Fr(k) 
jf(fc) = J drW s (k,r)j l (kr) . 



(21) 



In the present paper we consider two secondary effects 
that correlate with lensing deflections: integrated Sachs- 
Wolfe (ISW) effect and the Sunyaev-Zel'dovich (SZ) ef- 
fect. Since higher order effects such as the kinetic SZ 



T ISW (n) 



dr&(r, nr) 



(22) 



The weight function associated with the ISW effect is 
given by 



W iiiW (k) 



-3fL 



(23) 



which can be used to calculated the correlation between 
lensing pontential and the ISW effect through equa- 
tion 
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The SZ effect leads to an effective temperature fluctu- 
ation in the CMB given by the integrated pressure fluc- 
tuation along the line of sight: 



T sz (n) = g(x) 



dra(r) ^—jn e (r)T e (r) 



m e c* 



(24) 



where ut is the Thomson cross-section, n e is the elec- 
tron number density, r is the comoving distance, and 
g(x) — xcoth(a;/2) — 4 with x = hv/ksTcMS, is the spec- 
tral shape of SZ effect. At Rayleigh- Jeans (RJ) part of 
the CMB, g(x) — —2. For the rest of this paper, we 
assume observations in the Rayleigh- Jeans regime of the 
spectrum; an experiment such as Planck with sensitivity 
beyond the peak of the spectrum can separate out SZ 
contributions based on the spectral signature, g(x) [Q. 

For the correlation between lensing and SZ effect, we 
follow the halo model approach of (23) (see (24) for further 
details) which allows a semi-analytical approach to calcu- 
late the power spectrum of large scale structure pressure 
fluctuations. At RJ part of the frequency spectrum, the 
SZ weight function is 



W b/l {r) 



k B <7Tn e 
a(r) 2 m e c 2 



(25) 



where n e is the mean electron density today. With the 
halo model, we replace the clustering of dark matter with 
that of pressure when describing the SZ effect. The cross- 
correlation between lensing and SZ then involves the 
cross-power spectrum between pressure and dark mat- 
ter. 

In figure [l], we show the angular power spectrum 
of CMB anisotropies |25|] with secondary contributions 
through weak lensing, ISW and SZ effects. The SZ an- 
gular power spectrum was calculated using the halo ap- 
proach of [^H . The lensing angular deflection power spec- 
trum and the resulting correlation power spectra between 
lensing and ISW, and, lensing and SZ effects are shown 
in figure 0. 
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III. LENSING CONTRIBUTION TO CMB 

In order to derive weak lensing contributions to the 
CMB trispectrum, we follow Hu ||] and Zaldarriaga 
We formulate the contribution under a flat sky approxi- 
mation; this formulation is adeqaute given that we are 
mostly interested in non-Gaussian effects due to lens- 
ing at small angular scales corresponding to multipoles 
> 1000. In general, the flat-sky apporach simplifies the 
derivation and computation through the replacement of 
mode-coupling Wigner symbols through angles. 

As discussed in prior papers [§|||H > weak lensing maps 
temperature through the angular defelctions resulting 
along the photon path by 



e'(n) = 9(n + V<£) 

= e(n)+V^(n)V ,; e(n) 

+ iv ! 0(n)V J 0(n)V l V-'e(n) 



(26) 



As expected for lensing, note that the remaping conserves 
the surface brightness distribution of CMB. Here, O(n) 
is the unlensed primary component of CMB and 0*(n) is 
the total contribution. It should be understood that in 
the presence of low redshift contributions to CMB result- 
ing through large scale structure, the total contribution 
includes a secondary contribution which we denote by 
s (n). Since weak lensing deflection angles also trace 
the large scale structure at low redshifts, secondary ef- 
fects which are first order in density or potential fluc- 
tuations also correlate with the lensing deflection angle 

Taking the Fourier transform, as appropriate for a flat- 
sky, we write 



6(li) = / dhO(h)e- 



9(10 



e(ii')£(ii,ii'), 



(27) 



where 



i - 



d 2 h 



C 



(30) 



(2tt) 



c$_ hl ct;i(i-h)-hY 



As written, the second term shows the smoothing behav- 
ior of weak lensing through a convolution (eon. [30]) of 
the CMB power spectrum (see discussion in g). With 
respect to lensing contribution, there are two limiting 
cases: when 1 — li ~ 1 and cmb power is constant, one 
can take the CMB power spectrum out of the integral in 
the second term such that 



(27r) 2 C Ix ^ l > 



(31) 



Thus, there is a net cancellation of terms involving lens- 
ing potential power spectrum and Cf « Cf producing 
the well known result that lensing shifts but does not cre- 
ate power on large scales ||] . On small scales where there 
is no or little intrinsic power in the CMB, the second term 
behaves such that 



Cf 



(32) 



Here, the power generated is effectively the lensing of 
the temperature gradient associated with the damping 
tail of CMB anisotropy power spectrum. This small an- 
gular scale limit and its uses as a proble of large scale 
structure density power spectrum and mass distribution 
of collpascd halos such as galaxy clusters is considered in 



B. Trispectrum 

The calculation of the trispectrum follows similar to 
the power spectrum. Here, we explicitly show the calcu- 
tion for one term of the trispectrum and add all other 



terms through permutations. First we consider the cumu- 
x 0*(li" + 1/ — li) (li" • li')(li" + li' — li) • plants involving four temperature terms in Fourier space: 



We define the power spectrum and the trispectrum in 
the flat sky approximation following the usual way 



(^(lije^la)) = (27r) 2 Ml 12 )(7 ( 



i ' 



(e^io.-.e^u))^ 



(e*(ii) . . . e t (i 4 )) c = (27r) 2 <5 D (i 123 4)f (i!, i 2 , i 3 , h) . 



(29) 



A. Power spectrum 

The power spectrum, according to the present formu- 
lation, is discussed in M and we can write 



x 6(1 2 ) 



d% 



■Q(h')L(h,i 2 ') ) e(i 3 )e(u) 



7 ^ea/)^, o/A|e(i 2 W 2 ,i 2 ') 
xe(i 3 )e(i 4 )' 



(33) 



5 



As written, to the lowest order, we find that contributions 



come from the first order term in L given in equation (29) 
We further simplify to obtain 



(e'oo...^^))^ 



(27T)' 



epi'Mii - ii') (ii - ii') • ii' 



W) 2 



e(i 2 Xi 2 - 1 2 ') (i 2 - 1 2 ') • i 2 'e(i 3 )e(i 4 ) 



cZcf^h + 13)0(12 + h)) (h + 1 3 ) • 1.3 (i 2 + U) • i 4 



Perm. 



(34) 



where there is an additional term through a permutation 
involving the interchange of li + 13 with li + 14. Introduc- 
ing the power spectrum of lensing potentials, we further 
simplify to obtain the CMB trispectrum due to gravita- 
tional lensing: 



T e (li,l 2 ,l 3 ,l4 



q?f+i3i(ii + i3)-i 3 (ii + i 3 )-i4 



+Cg +hl (h + I 3 ) • ls(l 2 +I3) • 14 



+ Perm. , 



(35) 



where the permutations now contain 5 additional terms 
with the replacement of (I3, h) pair by other combination 
of pairs. 
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FIG. 4. The ratio of full variance, including non-Gaussian 
covariance, to that when non-Gaussian trispectrum contribu- 
tions are neglected. This ratio shows the fractional change in 
the variance, or the errors along the diagonal of the covariance 
marix. Out ot multipoles of a few thousand, the increase due 
to lensing is less than few percent and at the smallest scale 
there is no change from the lensing alone trispectrum as its 
contributions are substatially small. The lensing-SZ trispec- 
trum only makes noticable contributions at the smallest scale. 
The resulting changes from the lensing-ISW trispectrum is be- 
low 10~ 6 and can be ignored. 
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FIG. 3. CMB anisotropy trispectrum resulting from lens- 
ing, lensing-ISW and lensing-SZ correlations. The lensing 
trispectrum generally follows the shape of the CMB power 
spectrum while lensing-ISW and lensing-SZ trispectra depicts 
the behavior of correlation power between lensing and these 
secondary effects. 

The non-Gaussian contribution to the trispectrum, 
through coupling of lensing deflection angle to secondary 
effects, can be calculated with the replacement of 6(1,3) 
and 8(14) in equation ( |33"| ) by G s (l3) and O s (l4) contain- 
ing the sources of secondary fluctuations. Thus, we can 
no longer consider cumulants such as (9(li')9 s (l3)) as 



the secondary effects are decoupled from recombination 
where primary fluctuations are imprinted. However, con- 
tributions come from the correlation between s and the 
lensing deflection <fi. Here, contributions of equal impor- 
tance come from both the first and second order terms 
in L written in equation (^9|). First, we note 

(e(i 1 )...e(i 4 )) c = 

X (S^ 0(l2 ' )i(l2 ' l2,) ) 0S ( 1 3) 0S ( 1 4)) 

= -c h (L(i 2 ,-i 1 )e s (i 3 )e s (i 4 )) 
-c l2 ( J L(i 1 ,-i 2 )e s (i 3 )e s (i 4 )) 

+ / ^c ri (L(i 2 ,-i 1 ')L(i 1 ,i 1 ')e s (i 3 )e s (i 4 )) 

(36) 

Contributions to the trispectrum from the first two terms 
come through the second order term in L, with the two 
(j) terms coupling to s . In the last term, contributions 
come from the first order term of L similar to the lensing 
alone contribution to trispectrum. 

After some straightforward simplifications, we write 
the connected part of the trispectrum involving lensing- 
secondary coupling as 



G 



f e (l 1 ,l 2 ,l 3 ,l4) = 



cf s q 
13 1 



C«(l 3 -l 2 )(l 4 -l 2 ) 

f 1 4 ■ (ll + U)la ■ (1 2 



l3)Cg 



+ 1 3 • (ll 

+ Perm. . 

Note that the first two terms come from the first and 
second term in equation (36), while the last two terms 
in above are from the third term. As before, through 
permutations, there are five additional terms involving 
the pairings of (£3,^4). 

10° 



IV. POWER SPECTRUM COVARIANCE 

For the purpose of this calculation, we assume that 
CMB power spectrum will measure binned logarithmic 
band powers at several i^'s in multipole space with bins 
(37)of thickness 5L. 
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FIG. 5. The correlation coefficients of (a) lensing and (b) 
lensing-SZ contributions to the covariance of the CMB angu- 
lar power spectrum. In the case of lensing, the correlations 
shows the general behavior of lensing effectt on CMB where 
power is transferred from large angular scales with acoustic 
peaks to small angular scales, thereby correlating small and 
large angular scales at the tens of percent level or more. For 
experiments, such as Planck, that will measure the power 
specrum to multipoles of ~ 2000, the resulting correlations 
between bins are less than a percent. 



r d 2 i i 2 



(38) 



where A s (li) = J d 2 l is the area of 2D shell in multipole 
and can be written as A s (Zj) = 2irli8li + ir(5li) 2 . 
We can now write the signal covariance matrix as 



ij 



1 

A 



(2-Kf 



Aai 



2C 2 
d 2 l 



l 2 l 2 

3 l i l 3 



A sj {2-Kf 



T e (] — 1 1 — 1 



(39) 
(40) 



where A is the area of the survey in steradians. Again 
the first term is the Gaussian contribution to the sample 
variance and includes, in addition to the primary compo- 
nent, contribution through lensing and secondary effects. 
The second term is the non-Gaussian contribution. A 
realistic survey will also include instrumental noise con- 
tributions and we can modify the Gaussian variance to 
include the noise through an additional noise contribu- 
tion to the power spectrum 



cr 



(41) 



where JVj is the power spectrum of detector and other 
sources of noise introduced by the experiment. 

For the power spectrum covariance, we are interested 
in the case when I2 = — li with |li| = and I4 = — 13 with 
I X3 j = lj. This denotes parallelograms for the trispectrum 
configuration in multipolar or Fourier space. 

In the case of lensing contribution to the trispectrum, 
with the configuration required for the power spectrum 
covariance, we can write 



T e (V 
Cf.C^ 

+ cfcf 



2C ; e C ; e 



lli+ljl 

■4>4> 



[(l i + l j ).l i ] 2 + C | tf_ lj| [(l i -l J -)-l i ] 



\n-h\ 



-c 



(42) 



and includes all terms with no additional permutations. 

Similarly, for the lensing-secondary trispectrum we 
have 



T e (l»,- 
2 (L • 1, 



lj, lj, lj) — 
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2[l i -(l j -l i )] [^.(Ij-l^CfXty^ 

2 [i, •(!, + !,■)] [i r (u + i 3 )}cf;cf;c^ +1A 



(43) 



V. RESULTS 



In Fig we show the scaled trispectra, where 



A 

sq 



(0 



2tt 



T e (l,-I,Ix,-U 



,1/3 



(44) 



and l± = I and 1 • Ij_ = 0. In the case of lensing alone, 
the trispectrum is proportional to square of the CMB 
anisotropy power spectrum (see, equation |35| ) and the 
sharp reduction in power at multipoles greater than a 
few thousand is effectively due to the decrease in pri- 
mary anisotropy power at small angular scales. In the 
case of lensing-secondary correlation, the trispectrum is 
only propertional to one power of the CMB anisotropy 
power spectrum. Thus, the trispectrum now depicts 
more of the behavior of lensing-secondary correlation 
power spectrum shown in figure |[ The sharp decrease 
in the lensing-ISW trispectrum compared to that of the 
lensing-SZ effect is due to differences in the small angular 
scale power associated with lensing-ISW and lensing-SZ 
correlations. 

We can now use this trispectrum to study the contri- 
butions to the covariance. In Fig. [|, we show the ratio of 
the diagonal of the full covariance to the Gaussian vari- 
ance with the non-Gaussian term neglected: 



(45) 



and for bands U given in Table I. Here, we have used 
rather wide bins in multipolesuch that bin width is con- 
stant in logrithmic intervals in multipole space. This is 
the same binning scheme used by |2(| on 6° x 6° fields to 
investigate weak lensing covariance and later adopted by 
(□J. The two lines show the ratio when trispectra due 
to lensing and lensing-SZ correlations are used to calcu- 
late the covariance, resepectively. The square root of the 
ratio is roughly the fractional change induced in errors 
along the diagonal resulting from the non-Gaussian co- 
variance contributions. We do not show the ratio due 
to lensing-ISW trispectrum as the resulting changes are 
less than 10 -6 at all multipoles of interest. As shown in 
figure [|, the ratio is less than 20% for weak lensing and 
peaks at multipoles ~ 4000 while the ratio increases to 
smallest scale with the lensing-SZ trispectrum. 
The correlation between the bands is given by 



Cij = 



Ci 



yj Cii C 



(46) 



j j 



In Table I we tabulate the correlation coefficients for the 
CMB binned power spectrum measurements. The up- 
per triangle here is the correlations under the lensing 
trispectrum while the lower triangle shows the correla- 
tions found with the trispectrum due to lensing-SZ corre- 
lations. In the case of lensing contribution to the tripsec- 
trum, correlations depict the general shape of the CMB 
power spectrum while in the case of lensing-SZ contri- 
bution to the covariance, the correlation coefficients are 
more consistent with the shape of the lensing-SZ power 
spectrum. 

In Fig. ||, we show the correlation coefficients for (a) 
lensing and (b) lensing-SZ contributions to the covari- 
ance. Here we show the behavior of the correlation co- 
efficient between a fixed L (as noted in the figure) as a 
function of £j. Note that when = lj the coefficient is 
1 by definition; we have not included this point in the 
figure due to the apparent discontinuity it creates from 
the the dominant Gaussian contribution at U = lj . 

To better understand how the non-Gaussian contribu- 
tion scale with our assumptions, we consider the ratio 
of non-Gaussian variance to the Gaussian variance (see, 

00) 



= 1 + R, 



with 



R 



A T e 
(2tt)22C| 



(47) 



(48) 



In the case of lensing alone contribution to CMB trispec- 
trum, we can simplify the expression for R by noting 
that 



T e (l- -1 1 -1 



(49) 



where for an approximation we have taken lj • L = 0. 
Replacing the averaging of the product of {Cf.) 2 Cf^ with 
the product of two averages, we can simplify the ratio of 
TfJCf to obtain 



R ~ Ak5h 



itcf 



2tt 



(50) 



where < ... >Ai represents the averaging of the lensing 
potential powerspectrum, weighted by a factor of if to 



represent the deflection angles. The equation ([30|) repre- 
sents the general behavior of the non-Gaussian contribu- 
tion to the lensing trispectrum. The relative contribution 
from non-Gaussianities scale with several parameters: (a) 
increasing the bin size, through 5li (oc A s i), leads to an 
increase in the non-Gaussian contribution linearly while 
(b) the contribution is determined by the shape of the 
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lensing potential power spectrum. At large multipoles, 
when Cf^ is a rapidly decreasing function, there is no 
significant contribution to the trispectrum and thus to 
the covariance. This is contrary to the general assump- 
tion that lensing is a small scale phenomena and that 
small-scale signal in CMB should be strongly affected by 
the non-linear nature of weak lensing. In the case of 
CMB, however, lensing effect is a large scale phenomena, 
especially given the power spectrum of lensing deflection 
angles peak at multipoles of ~ 40 and most of the lensing 
power is concentrated at degree angular scales instead of 
arcsecond scales. 

For upcoming wide- field experiments, especially those 
involving satellite missions such as MAP and Planck, we 
do not expect non-Gaussianities to limit the interpreta- 
tion and the cosmological parameter extraction from the 
measured CMB power spectrum. For these wide-field ex- 
periments, the width of the bin in multipole space will 
be of order at most few tens; for such small bin widths in 
multipole space at I ~ few hundred will lead to a signifi- 
cantly reduced ratio of non-Gaussian to Gaussian contri- 
bution from what we have considered where SI ~ I. Also, 
note that the cosmologically interesting acoustic peak 
structure and the damping tail of the CMB anisotropics 
is limited to multipoles below ~ 1000. In this range, 
there is no significant non-Gaussian contribution related 
to lensing and lensing-secondary correlations. There are, 
however, ground-based experiments (e.g., Cosmic Back- 
ground Imager [E8|) for which the non-Gaussianities due 
to lensing may be important. These small angular scale 
experiments, which probe the anisotropy power between 
multipoles of ~ 1000 and 4000 or so are likely to be lim- 
ited to small areas on the sky and will utilize wide bins 
in multipole space when estimating the power spectrum 
in order to increase the signal-to-noise associated with 
its measurement. In such a scenario, it may be necessary 
to fully account for the full covariance when interpreting 
the power spectrum at small angular scales. 

VI. SUMMARY 

The upcoming small angular scale CMB anisotropy 
experiments are expected to provide first measurements 
of the power spectrum related to the damping tail and 
secondary anisotropics. At such scales, important non- 
linear effects and secondary contributions can imprint 
non-Gaussian signals on the CMB temperature fluctu- 
ations. Here, we discussed one aspect related to the 
presence of non-Gaussianities when measuring the CMB 
anisotropy power spectrum involving a contribution to 
the covariance resulting from the higher order four-point 
correlation function, or a trispectrum in Fourier space. 
Here, we discussed the non-Gaussian contribution to the 
power spectrum covariance of CMB anisotropies resulting 
through weak gravitational lensing angular deflections 
and the correlation of deflections with the integrated 



Sachs- Wolfe effect and the Sunyave-Zel'dovich effect. 

With substantially wide bins in multipole space, the 
resulting non-Gaussian contribution from lensing to the 
binned power spectrum variance is insignificant out to 
multipoles of few thousand containing acoustic peaks and 
the damping tail, which are of substantial interest for cos- 
mological parameter estimation purposes. For upcoming 
satellite based near all-sky experiments, we do not ex- 
pect non-Gaussianities to limit the cosmological parame- 
ter extraction from CMB power spectrum measurements 
and their interpretation. For small angular scale ground- 
based experiments with substantially limited sky cover- 
age, however, the presence of non-Gaussianities should be 
accounted when interpreting any measurements at angu- 
lar scales corresponding to fewe arcminutes or multipoles 
~ 4000. Observational results related to this part of the 
CMB power spectrum is expected from a wide-variety 
of experiments based on interferometric and bolometric 
techniques. 
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